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Abstract 

We consider the large- N $ 4 theory with spontaneously broken symmetry at finite 
temperature. We study, in the large-N limit, quantum states which are character- 
ized by a time dependent, spatially homogeneous expectation value of one of the 
field components, 0jv(i), and by quantum fluctuations of the other N — 1 com- 
ponents, that evolve in the background of the classical field. Investigating such 
systems out of equilibrium has recently been shown to display several interesting 
features. We extend here this type of investigations to finite temperature systems. 
Essentially the novel features observed at T = carry over to finite temperature. 
This is not unexpected, as the main mechanisms that determine the late-time 
behavior remain the same. We extend two empirical - presumably exact - rela- 
tions for the late-time behavior to finite temperature and use them to define the 
boundaries between the region of different asymptotic regimes. This results in a 
phase diagram with the temperature and the initial value of the classical field as 
parameters, the phases being characterized by spontaneous symmetry breaking 
resp. symmetry restoration. The time evolution is computed numerically and 
agrees very well with the expectations. 
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1 Introduction 



The investigation of the O(N) vector model at large N has a long-standing history 
in quantum field theory |2|, [|. One of the main aspects was the question of 
symmetry restoration at high temperature that for some time was controversial. 
The dynamical exploration of a special class of nonequilibrium properties has 
been developed only recently ||, |,|(| [/] . 

The out-of equilibrium configuration that has been studied mainly is char- 
acterized by an initial state in which one of the components has a spatially ho- 
mogeneous classical expectation value <f)(t). This implies that the other N — 1 
components ^(x, t), i = 1 . . . N — 1 have a mass that is different from the mass in 
the ground state. This means that their initial state is related to the Fock space 
vacuum state by a Bogoliubov transformation. The evolution of the system is 
governed by the classical equation of motion for the field and by the mode 
equations for the quantum fields ip(x,t). The expectation value (^(x, t)^(x, t)) 
appears in both equations of motion, this constitutes the quantum back reaction. 
In the one-loop approximation, in contrast to the large-iV approximation, this 
quantum back reaction only appears in the classical equation of motion. This 
leads to decisive differences in the late time behavior. 

We have previously || carried out such dynamical computations for the O(N) 
vector model in the limit of large N at finite temperature for the case of unbroken 
symmetry, i.e., with a positive mass term. Here we will consider the case of 
spontaneously broken symmetry. In this case, at low temperatures the fields 
ipi(x,t) will be the Goldstone modes. This is the case for the ground state at 
T = and at finite temperature; for nonequilibrium initial states these modes 
become massless when the system settles to a stationary state at late times. 
Symmetry restoration happens at high temperature and at large values of the 
initial field 0(0); then at late times these modes stay massive while the classical 
field vanishes, and thereby the spontaneous symmetry breaking disappears. 

Our investigation, as well as the analogous ones at T = 0, are limited to fields, 
masses (as solutions of the gap equation) and temperatures much smaller than 
the scale of the Landau ghost m x = mi exp(87r 2 /A), where mi is a renormaliza- 
tion scale, taken of order y/\v. So the question of symmetry non-restoration at 
"really" high temperatures || will not be addressed here. 

The plan of the paper is as follows: in section 2 we introduce the model 
and set up the equations governing the nonequilibrium evolution. In section 3 
we discuss the renormalization of the equations of motion and of the energy- 
momentum tensor, some details are referred to Appendix A. In section 4 we 
discuss the phase structure of the system as a function of temperature and initial 
conditions. In section 5 we present the results of the numerical computations. 
Some conclusions are drawn in section 6. 
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2 Formulation of the model 



We consider the O(N) vector model with the Lagrangian 



(2.1) 



where = 1,..,N are N real scalar fields. The nonequilibrium state of the 
system is characterized by a classical expectation value which we take in the 
direction of (p^. We split the field into its expectation value (f) and the quantum 
fluctuations ip via 

0*(x, t) = SifVNW) + ^(x, t) . (2.2) 

In the large- N limit one neglects, in the Lagrangian, all terms which are not of 
order N. In particular terms containing the fluctuation ip^ of the component 
4>n are at most of order VN and are dropped, therefore. The fluctuations of 
the other components are identical, their summation produces factors N — 1 = 
N(l + 0(1/N)). In the broken symmetry case these are the Goldstone modes. 
Identifying all the fields ipi, ..ipN-i as ip the leading order term in the Lagrangian 
then takes the form 

C = N (£ + + d) , (2.3) 



with 



£<t> 



(2.4) 
(2.5) 
(2.6) 



where ip 2 is to be identified with J^V'VV-^- 

We decompose the fluctuating field into momentum eigenfunctions via 



d 3 k 



(2vr) 3 2^ 



a k U k (t)e^ + a{U* k (t)e 



-ikx 



(2.7) 



with Uko — y m o + k 2 . The mass m will be specified below. This field decompo- 
sition defines a vacuum state as being annihilated by the operators a^. 

The equations of motion for the field <p(t) and of the fluctuations Uk{t) have 
been derived in this formalism by various authors || [K], [U] . 



We include in the following the counterterms that we will need later in order to 
write the renormalized equations. The equation of motion for the field <fi becomes 



<j>(t) + 8m 2 (j>(t) - \v 2 <p(t) + (A + 5\)(j>(t) (j> 2 (t) + T(t, T) 







(2i 
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Here J-'it, T) is the divergent fluctuation integral; it is given by the average of 
the fluctuation fields defined by the initial density matrix. For a thermal initial 



state of quanta with energy uj k0 = Jk 2 + m 2 , it is given by 



d 3 k 



(2tt)32^ 

The mode functions satisfy the equation: 

d 2 



coth 



U k {t)\< 



dt 



U k (t) = 



(2.9) 



(2.10) 



and the initial conditions 

U k (0) = l ; U k (0) = -iu ko . 
The time dependent frequency ui k {t) is given by 

uj 2 (t) = k 2 + M 2 (t) 

with the time dependent mass 

M 2 {t) = -Xv 2 + 5m 2 + (A + 5X) [(f) 2 (t) + F(tj . 

Using this definition the classical equation of motion can be rewritten as 

0(t)+M 2 (t)(f)(t) =0 (2.14) 

which is the same equation as the one for Uk(t) with k = (zero mode). Of 
course the initial conditions are different and <f>(t) is real. 

As in our previous work we rewrite the mode equation in the form 



(2.11) 



(2.12) 



(2.13) 



d 2 2 

dT 2+UJk0 



U k (t) = -V{t)U k {t) 



(2.15) 



whereby we have defined the time-dependent potential V(t) = M. 2 (t) — .M 2 (0); 
we further identify m = M. (0) as the "initial mass" . 

The average of energy with respect to the initial density matrix is given by 



^(f)^(-y + 5mV(t) + - 



A ~ OA ,A, s - . 

(f)\t) + 5 A 



+ 



d 3 k (3u k0 f 1 

- — — coth < - 

{27c) 3 2u k0 2 1 2 



U k (t)\ 2 + -u 2 k (t)\U k (t)\ 2 



(2.16) 



A + 5X _ 2 



4 



T 2 {t,T) . 



3 Note that twice the last term, with positive sign, is included in the fluctuation energy, since 
\(t) contains F(t,T). 
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It is easy to check, using the equations of motion ( [2.14j ) and ( |2.1Q| ) , that the energy 
is conserved. The energy density is the 00 component of the energy-momentum 
tensor. The average of the energy momentum tensor for our system is diagonal, 
its space-space components define the pressure which is given by 



p = <f) 2 {t)-£ + 5t 
d 3 k 



d 2 



+ 



/ 



dt 2 
coth 



<p 2 (t)+F(t,T) 



+ 3-) \Uk(t)\' 



(2.17) 



5£ is the renormalization of the conformal coupling term £ (g^d 2 — d^d u )<j) 2 , which 
has been used for the improved energy momentum tensor fT2"]. 



3 The renormalized equation of motion 

The expressions for the time-dependent mass A4 2 (t), the energy density £{t) and 
the pressure are still undefined as they involve divergent integrals over the fluctu- 
ations. Our approach to regularization and renormalization has been presented 
previously |13|, ||. It is based on expanding the fluctuations Uk(t) and sub- 
sequently the various integrals involving these fluctuations with respect to the 
time-dependent potential V(t). As this procedure has been presented elsewhere 
in detail we just give the outline, here. 

The expansion of the fluctuations with respect to V(t) is given in Appendix 
A. We use this perturbative expansion in order to single out the divergent con- 
tributions in the fluctuation integral. One finds 



F(t) = J_i(mo,T) - i_ 3 (m ,T) [M 2 (t) - M 2 (0)} +F Rn (t,T) 
where the finite part of T{t, T) can be written as 



(3.1) 



d 3 k 



+ 



(2vr)3 4u 3 
d 3 k 1 



dt' cos [2oo ko (t - t')] V(t') coth 



2Re/f + 



(2tt)3 2u k0 

and where the divergent integrals are defined as 



coth 



I3u k0 



J_i(m ,T) 



d 3 k 



(2tt) 3 2oo ko 
d 3 k 1 
(2vr) 3 4a; 3 



1 + 
1 + 



2 



(3.2) 



1 



I-i{m ) + S-i(m , T) ,(3.3) 
/_3(m ) + S_ 3 (m ,r).(3.4) 
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The integrals J_fe(mo) are the genuine divergences which appear in the renormal- 
ization at T = 0. Their dimensionally regularized form is given by 




-mS/_ 3 (mo) - ^5 • (3.6) 

The additional temperature dependent terms £_fc(m ,T) are finite. They are 
defined as 

f d 3 k 1 

S - l(m °' T) = J {2,f ^ (e^o _ 1) > ( 3 - ? ) 

S - 3(m °' T) = / (2tt)3 2^(^-1) - (3 - 8) 

It is convenient to include these finite terms into the definition of J r fi n (t, T). Then 
the time dependent mass takes the form 

M 2 (t) = X(<p 2 -v 2 )+5X<p 2 +5m 2 + (X+5X) [l-i(wo) -/_ 3 (w )V(t) +^an(t,T)] , 

(3.9) 

with 

J- fin (t,T) = S_ 1 (mo,r)-V(t)S_ 3 (mo,T) + ^ fin (t,T) . (3.10) 

The time dependent mass ( |3.9| ) contains both renormalization constants 5m and 
<5A. Furthermore, its definition by this equation is implicit, Ai 2 (t) also appears 
on the right hand side of ( |3.9| ) in V(t). 

We now have to fix the renormalization counterterms in such a way that the 
relation between the time dependent mass and 4>(t) becomes finite. An additional 
constraint derives from the requirement that the renormalization counterterms 
should not depend on the initial condition, but only on the parameters appearing 
in the Lagrangian, i.e., A and v and renormalization conventions. 

We first determine SX by considering the difference 

V(t) = M 2 {t)-M 2 {0) (3.11) 
= (X + 6X) [0 2 (t)-0 2 (O)-/_ 3 (m o )V(t)+ ^ fin (t,T)-J- fin (0,T)" . 

The divergent parts depend on the initial mass mo. We have to replace this 
by a renormalization scale independent of the initial conditions. In Ref. || 
we had chosen the scale m, where m was the mass parameter appearing in the 
Lagrangian. Here the analogous mass squared would be m 2 = —Xv 2 and so m 
would be imaginary. We therefore choose another scale m\ which we do not 
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specify here. In the numerical computations we have used the physical Higgs 



mass m 



m 



H 



2\v 2 . 



We rewrite the implicit equation for V(t) as 
V(t) [1 + (A + 5A)7-3K)] = (A + SX) [<P 2 (t) - 2 (O) - [/_ 3 (m ) - /_ 3 (mi)] V(t) 



and require 



+^fi„(t,T) - Jfi n (0,T)} 

A + 5A 



1 + (A + 5A)/_ 3 (m x ) 
Solving with respect to 5X we find 

A 2 J- 3 (mi) 



A . 



5\ 



1 - A/_ 3 (mi) ' 
Inserting this relation into ( ft . 1 2| ) we find 

V(t) = AC [0 2 (t) - 2 (O) + F &n (t, T) - J- fin (0, T) 

with 

C 



1 



1 



1 + A [/_ 3 (m ) - J_ 3 (mi)] 









H 


- A > 




16tt 2 


Wo/ 



(3.12) 
(3.13) 

(3.14) 

(3.15) 
(3.16) 



Eq. (|3.15| ) is a finite relation for the potential V(t) since the difference [J_3(mo) — 
J_ 3 (mi)] is finite. Going back to Eq. ( ft.lOQ we realize that jF fin on the right hand 
side contains itself a term proportional to V(t). Taking account of this term we 
rewrite V(t) in terms of J 7 ^ as 



with 



C 



V(t) = XC T [<p 2 (t) - </> 2 (0) + Fsn(t,T) 
1 



A 



16vr 2 



In 



m\ 



(3.17) 
(3.18) 



+ A£_ 3 (m ,T) 



Recall that .Ffin(£) is the mode integral of second order in V(t) and vanishes at 
t = 0. 

We now go back to equation ( |3.9| ) which we take at the initial time t — 0: 

m 2 = M 2 (0) = A[0 2 (O) - v 2 } + 5X(f> 2 (0) + 5m 2 + (A + 5X) [l_i(m Q ) + ^(0, T) . 

(3.19) 

This is an implicit relation between mo and 0(0) which, however, contains still 
the infinite quantities SX, 5m and J_i(m ). Using Eq. (|3l|) we can rewrite Eq. 
( KW) as 



m r 



Xv 2 + 6m 2 ) + (A + SX) 



2 (O) - m 2 J_ 3 (m ) - 



167T 



+ ^fin(0,T) 



(3.20) 
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As renormalization condition we require ttlq to vanish, for temperature T = 0, 
at the minimum of the potential <f> — v, as it is the case on the tree level. We 
note that m 2 , = is not the curvature of the tree level potential at <ft = v 
which is m 2 H = 2Xv 2 . It is the mass of the fluctuations at <f> = v in the large- N 
approximation. For T = we have ^ n (t — 0, T — 0) = T l _i(m ,T = 0) = 0. 



Setting m = 0, 0(0) = v in the gap equation ( |3.20| ) we get immediately 



5m 2 



-5\v 2 



A 2 f 2 /_ 3 (mi) 
I-AL3N 



(3.21) 



Inserting this into Eq. ( |3.20| ) we obtain the renormalized gap equation 



777.Q = AC 



(0)-t; 2 -— iL + S _ l(m0;T ) 



16tt 



(3.22) 



For the numerical computation it is easier to choose some m 2 , > and to use the 
gap equation solved for 2 (O): 



\0) = ^ + v 2 + ^- 2 
A 16n 2 



fl+ln^| ] -E_ a (mo,T) . (3.23) 



For t > the renormalized relation for the mass squared M. 2 (t) we find, using 
Eqns. ([TT5D and (|3722|) , is 



>J 2 (t) = m 2 + V(t) = AC 



m, 



167T 



^ + J- fin (t,T) 



(3.24) 



Having thus obtained a finite relation between <f>(t) and the equations of 

motion for the classical field <f>(t) and for the modes Uk(t) are well-defined and 
finite. 

The way in which we have renormalized has made the cutoff disappear. This 
was possible only to the extent that we could safely neglect corrections of order e 
in the evaluation of the divergent integrals. One way of achieving this is to take 
the limit e — > 0. This implies for the bare coupling Aq 



Ao 



lim 

e^0 



A 



A 2 



(T 



(3.25) 



IQtt 2 e 



so this is the case of "negative bare coupling" as discussed in [0. One can leave 
the cutoff finite, however, as long as the masses and momenta are much smaller 
than the scale of the Landau ghost, m x = m\ exp(87r 2 /A). This will be case 
here. This is not related to a pragmatic momentum cutoff that we apply to the 
convergent integrals of the finite part. 
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While we have found here the gap equation as a self-consistency condition, it 
can also be derived |5], [7j from a potential (free energy) which here takes the form 



V{ml®\T) 



+ 



rrir 



— v — 



+ 



rrir 



2X 32vr 2 



In 



m5 



d 3 k 1 



In [1 - exp(-/3w )] 



(2tt) 3 /5 

The gap equation then follows from the condition 

dV(ml,(f) 2 ,T) 



(3.26) 



dmk 







(3.27) 



It should be mentioned here that the gap equation has two solutions, one of which 
lies above the scale of the Landau ghost, m x = mi exp(8ir 2 / X 2 ) . In the sense that 
we consider here the model as giving rise to a low energy effective theory we 
discard this high mass solution, and its discussion. The solution we consider is 
the low energy one which is of order y/Xv or m±. 
The energy density is given by 



£ = h 2 (t) + \(X + 5X)(<p 2 -v 2 ) 2 + 5A 



-£a(t,T) 



X + 5X 



T 2 {t,T) . 



(3.28) 



Here we have used already that 5m 2 = —5Xv 2 , and part of the "cosmological 
constant" counterterm 8 A is included in SXv 4 /4. The fluctuation energy is given 
by 



£a{t,T) 



d 3 k 



coth 



~\U k {t)\ 2 + -u>t{t)\U h {t)[ 



(27t) 3 2^ ■ 2 
We again split off the temperature-dependent contribution via 

£ & (t,T)=£ 6 (t,0) + A£ Q (t,T) , 
where the second term on the right hand side 



A£ fl (t,T) 



d 3 k 



(27r) 3 2cj fc0 e^w - 1 1 2 



l\Uk(t)\ 2 + l^t(t)\u k (t)\ 2 



(3.29) 



(3.30) 



(3.31) 



is finite. The divergences of the first term are given [U| by the decomposition 
ffl(t, 0) = h(m ) + ^V(t)/_i(m ) - iv 2 (t)/- 3 (mo) + ^W*, 0) (3.32) 
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with 



£a,fm(£, 0) — - 



d 3 k 



(2tt)32^ 



+ 



V(t) 



2Re/i 1} + \ti 



(1)|2 



+ 



V 2 (t) 



A:0 



(3.33) 

We denote the sum of £n t sn{t, 0) and A£ fl (t, T) finite contributions as £a_fi n (t, T). 
The expression for the energy then takes the form 



£ 



1-2 A + 5X , , 
V 4 



2 + £ fl , fin (t,T) +/i(m ) + -V(t)/_i(m ) - iv 2 (t)/_ 3 (mo) 



A + <5A _ 2 



^ 2 (t,T) + <5A . 



(3.34) 



In addition to the divergences arising from £a(t, T) we have to take into consider- 
ation those of J-" 2 (t,T) which we have analyzed above. If all divergences and the 
renormalization constant 5X are inserted, the expression turns out to be finite, 
i.e., the remaining counterterm 5 A is needed only for a finite renormalization. 
We require the energy to vanish at T = for <fi(t) = v, which implies mo = 0. 
Then 5A = 0. There remains a finite constant dependent on the initial condition 



AA 



1287T 2 



1 + 



and the energy is given by 

1 • 



£ 



A 



,2\2 



2 AC 

16tt 2 



-Am 2 



(3.35) 



(3.36) 



A 



+£ &M (t, T) - -CT\Jt, T) + AA 



Here Am 2 is given by 

Am 2 

We write the pressure in the form 



-AC 



16tt 2 



P = <P 2 (t)-£+ P{l (t,T) + 6£ 



dt 2 



\t)+F(t,T) 



(3.37) 



(3.38) 



The renormalization does not differ form the case of unbroken symmetry discussed 
in Ref. M and is not presented again. We find 



Aa; 



A/_ 3 (m) 



6(1 -As) 6(l-A/_ 3 (m)) 
The final result for the renormalized pressure reads 



p = <P 2 (t)-£+ PaM (t,T) 



rrir 



96vr 2 48tt 



96vr 2 



ml* 



ill 



(3.39) 



V(t) (3.40) 
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with 



Pfl,fin(*, 0) 




/ 1 m 2 \ 

+ T^X + orV cos(2^ fc0 t)V(0) (3.41) 




+ \fi 1 \t)\ 2 - 2Re [iwrfW + iWiPW^)] } • (3.42) 



4 Analysis of the gap equation and of the phase 
structure 

The dynamical evolution of the nonequilibrium system depends on two parame- 
ters, the temperature T and the initial amplitude of the classical field 0(0) = <fto 
which in analogy with thermal equilibrium systems can be considered as an ex- 
ternal parameter. There are two regions from which we can start the system, 
which we will call regions II and III. There is, in addition, one region into which 
the system can evolve when one considers <p(t = oo) = 0oo and not 0o as the 
external parameter. We call it region I. In this section we will characterize these 
regions and describe the dynamical evolution as to be expected from the analysis 
at T = 0. This analysis is based on certain empirical results || [7) that, though 
unproven, seem to be at least almost exact. We will generalize these results to 
finite temperature in a plausible heuristic way, to be confirmed by the numerical 
computations. We think that the way in which we generalize these results will 
give a further clue to understanding them. 

4.1 Region I, ml < 

The gap equation requires Wq to be positive. The point where Wq = marks an 
initial condition that leads to a solution = const., if 0o — as we will assume 
in the following. For T = this stationary amplitude is = t>. For T > we can 
easily find this amplitude as well. Indeed for m = the integral E_i(m ,T) is 
given by its value for massless quanta, i.e., 



£-i(0,T) 



12 ' 



(4.1) 



therefore 



= 0oLo=O 




(4.2) 
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For 0o < 01 (T) the gap equation has no real solution tuq. The region below the 
boundary ( |4.2|) is region I. 

If nevertheless one wants to start the system with 0o in region I one faces 
the problem that in this region the gap equation requires to be negative. 
Then the low-momentum modes with k 2 < —ml have imaginary frequencies. So 
from an orthodox point of view (to which we adhere here) the system cannot 
be quantized properly. One may avoid this problem by redefining the dispersion 
relation for the initial frequencies via u>l Q = k 2 + \ml\ in this region. Of course 
at t > M. 2 {t) will be negative so the "mass squared" changes sign at T = 0, a 
situation called "quench" in analogy by a similar transition form a stable to an 
unstable state by a sudden drop of temperature or inversion of a magnet field. 
On the other hand the amplitude <p(t) can reach this region at late times, but 
then it is in a quantum state different form the ones we use as initial states. 

4.2 Region II: ml > 0, = 

We now assume O is started above the boundary value (|4.2|). If O is not too 
large the system may, at t > 0, enter a region where M. 2 (t) < 0, i.e., region 
I. Then the quantum fluctuations with momenta k 2 < —A4 2 (t) will increase 
exponentially, signalling instability. This causes M. 2 {t) to increase so that it is 
driven back to a value A4 2 (t) > 0. If the initial amplitude 0o is sufficiently 
small this forth- and-back reaction will lead Ai 2 (t) to stabilize at M. 2 ^ = 0. So 
at late times ^ n (t, T) is determined by quantum modes U^it) that oscillate with 
time-independent frequencies ujoo = k, it becomes stationary as well and will be 
positive. Therefore <fi(t) stabilizes at some value 

(ft = v 2 - .F fin (oo, T) < v 2 - ,F fin (0, T) . (4.3) 

This is entirely analogous to the behavior found at T = [0 . We call the region 
of initial values 0o leading to this late-time behavior region II. 

The stabilization by back-reaction onto the fluctuations obtained in the large- 
N approximation is not present in the one-loop approximation. In this approxi- 
mation, once <p(t) dips into the unstable region 4>{t) < v/V3, the mass squared 
of the fluctuations becomes negative and the low momentum modes evolve ex- 
ponentially. The effective mass of the classical field increases exponentially as 
well, and continues to do so, but the mass squared of the fluctuations stays nega- 
tive. The amplitude <p(t) is driven towards zero. Nevertheless the classical energy 
continues to increase, as the field oscillates faster and faster, this energy being 
extracted from the energy of the quantum fluctuations. Obviously this signals 
the instability of the quantum vacuum, as already apparent from the fact that the 
effective potential is complex in this region. We will illustrate this by a numerical 
example, to be presented in the next section. 
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At T = the final value 0oo was found to be related to the initial value <fio by 
an empirical relation 



%(2v 2 -cf> 2 ) T = 0. (4.4) 



It is not obvious how to generalize this relation to finite temperature. It was 
remarked in Ref. ]/] that the relation only depends on the initial, purely classical, 
energy, which is given by E = A(0 2 — t> 2 ) 2 /4. Obviously it satisfies the constraints 
that 0^ = v 2 if 0o = v 2 , and that (p^ = if classically the system can reach 
the maximum of the potential; this happens at 0q = </>|(T = 0) = 2v 2 . So Eq. 



(4.4) seems to be related to energy considerations. We further observe that the 

I = 2v 2 -<P 2 



classical turning point is at 4>q = 2v 2 — §\ so that one may write Eq. (|4.4| ) as the 



geometric mean 



& = vW§ • ( 4 - 5 ) 



This form turns indeed out to lead to the correct generalization for finite tem- 
perature. 

Obviously the relation is characterized by the motion at early times when the 
quantum fluctuations have not yet evolved. When discussing renormalization we 
have made an expansion with respect to the "potential" V(t) which vanishes at 
t = 0. So the same expansion can be used to study the early time behavior. In 
the energy the coefficients of the terms of first and second order in V have been 
absorbed into renormalization constants. However, the thermal fluctuations are 
not absorbed in this way and will add to the classical terms in an early time 
expansion. These appear in the energy, see Eq. (|3.34j) , via 

A£ fl (t,T) = Si(mo) + iv(*)£_i(mo) - ^V 2 (t)S_ 3 (m ) + 0(V 3 ) (4.6) 



as a part of £^ n {t,T) and via Eq. (|3.10| ) in jF fin (t, T). Taking these expansions 



4 



into account the energy can be written in the form 

a0 4 + a0o + H 2 + Hi + c0 2 0o] + const - ( 4 - 7 ) 
up to terms of order V 3 . We need the coefficients 

a = l-AC T S_ 3 (m ,T) , (4.8) 



b = -2 



f 2 -S-i(m 0) T) , (4.9) 



c = ACC T S_ 3 (m ,T) . (4.10) 
The classical turning point is given by 



- 2 _ b+(a + c)cj> 2 _ 1 



'0 



1 - AC T S. 



2v 2 - E_x - (1 + AC t E_ 3 )0q , (4.11 
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so that we are led to suppose 



0oo(T) = V l-AC T S 3 V^ ^ " 2S - X " (1 + XC ^)^] • (4-12) 

We find indeed (see below) that this relation is very well fulfilled numerically. 
According to this formula the region II is limited by the requirement that the 
expression in the square root be positive, so the boundary between region II and 
the new region III is given by 



, 2 = o v 2 -X-i(m ,T) 
2 l + AC T S_ 3 (m ,T) 



(4.13) 



We note that the relation is implicit, the value of mo that appears on the right 
hand side is related to 4>\ on the left hand side by the gap equation. 



4.3 Region III, = and M 2 ^ > 

If the value </>o becomes larger than 02 the stationary state with constant and 
vanishing mass A4 2 (t) is no longer attained, and the system reaches another 
asymptotic regime where Ai 2 (oo) ^ whereas <p(t) — ► 0. This regime is similar 
to the one that describes the late time behavior for the unbroken symmetry case. 
We call the region of initial values 0o that leads to such a behavior region III. 

There are two phenomena that characterize the transition to this region. On 
the one hand the stabilization of the system is taken over by the phenomenon 
of parametric resonance. On the other hand the system has enough energy so 
that 4>(t) can move over the maximum of the potential at = 0, and indeed will 
oscillate around = 0. Accordingly the threshold value of 0o at which these two 
phenomena set in can be characterized by two - a priori unrelated - criteria. Both 
rely on plausible assumptions, which at T = lead to the same prediction for 
the critical value of <pQ. 

The criterion based on the energy consideration has been presented in the 
previous subsection, we now describe the criterion supplied by the phenomenon 
of parametric resonance. For the case of unbroken symmetry it was found at 
zero |J and finite temperature |J, that the late time behavior is described by 
an empirical sum rule which relates Ai 2 ^ to the initial amplitude. For T = 
an analogous sum rule was found to hold for the case of spontaneously broken 
symmetry as well fTj. It is given by 

A4 = -1 + | ■ (4-14) 

Here \i and r\ are normalized in such a way that the classical equation of motion 
at early times, i.e., in the parametric resonance regime without back reaction, 
reads 

T]" -ri + r] 3 = , (4.15) 
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where the prime denotes a derivative with respect to r = at and where r\ = (3(f), 
also \i = M./a. With r](r) a solution of Eq. ( }4.15|) the mode equation becomes a 



Lame equation. The sum rule implies ||, that the frequencies uj 2 (t) = Ai 2 (t) + k 2 
are shifted outside the parametric resonance band of the Lame equation. Though 
there is no rigorous derivation for the sum rule, it accordingly seems related to 
the parametric resonance phenomenon. 

As the shift of the frequencies outside the parametric resonance region must 
have happened at the end of the phase where the evolution of the system is 
described by parametric resonance, we will again consider the initial classical 
evolution. Again, in addition to the classical terms we have to take into account 
the terms due to the thermal fluctuations. In terms of the parameters introduced 
in the previous section the equation of motion is given by 

'<p + XCcuf) 3 + ^C(b + c(f)l)(f) = . (4.16) 

Comparing to the normalized equation ( |4.15| ) we determine the factors a and (3 
to be 

a = + c4>\ , (4.17) 

2a 
b + c<pl 



P = \/-rr^s. ( 4 - 18 ) 



so that the asymptotic mass is given by 

Ml = a\-l+ l -^<Pl) (4.19) 

= Ac{^ 2 + E_i(mo,r) + i[l + ACrE_3(mo,T)]^} . 

Again <f) Q and m are related by the gap equation. At the transition from region 
II to region III the asymptotic mass vanishes. It is easily seen that this criterion 
leads to an identical equation for the boundary, i.e., Eq. ( |4.13| ). 

The field amplitude decreases to zero at late times, in this regime. So the 
symmetry is restored dynamically at high excitation characterized by a high 
value of 0o- 

At the critical temperature T = y/l2v both boundaries 4>\{T) and 4>2{T) 
become zero. Above Tc the behavior of the system is the same as for region III, 
for all initial values of 0q. While at the border between region I and II there was 
a lowest value for O f° r obtaining real solutions of the gap equation, now there 
is a lowest value of mo, the one for which 0o — 0. It is obtained by solving the 
gap equation for 0o — and agrees with the thermodynamical equilibrium value 
mp at that temperature, as defined, e.g., in Eq. (3.38) of Ref. 0. Of course with 
4>o = the system remains static. 
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Having denned the three regions by the two boundaries ( f4.2| ) and ( [4.13| ) we 
present, in Fig. 1, a phase diagram in the 0q — T plane. Fig. 2 shows the phase 
diagram in the ml — T plane, displaying, above Tc, the region m < mp which 
is excluded as an initial condition. We have to stress that the boundary between 
regions II and III relies on an empirical relation. 

The symmetry restoration above a critical temperature is expected naively. 
However, if the temperature becomes nonperturbatively large, T ~ \fYlm\ exp(87r 2 / A), 
the gap equation does not have solutions any longer. Then the free energy attains 
its maximum at the boundary mo = and the O(N) symmetry is again broken 
[[|. This phenomenon of "symmetry non-restoration", as well as the existence of 
the second solution of the gap equation above m x = m± exp(8ir 2 / X) will not be 
discussed here, as it is not part of the low energy effective theory. 



5 Numerical Results 

We have discussed already in the previous section the type of nonequilibrium 
behavior to be expected in the different regions of phase space. The numerical 
results follow these expectations. We have chosen generally the parameters v = 1 
and A = 1. We present results for the various regions in the T, <fio plane. The 
critical temperature is 2y / 3 = 3.464. We choose the temperatures between T = 1 
and 4, the latter one being above the phase transition. The numerical method has 
been described in M. We just recall that all the integrals computed numerically 
are finite, so cutting off the momentum integration at some reasonable value is 
unrelated to cutoffs used for renormalization. 

We first consider initial conditions in region II. The expectation value of 0, 
shown in Fig. 3 becomes constant and different from zero as t — > oo. This signals 
spontaneous breakdown of the O(N) symmetry. As displayed in Fig. 4 the 
mass Ai 2 (t) vanishes as t — > oo, as expected form the Goldstone theorem. The 
momentum distribution of the quantum fluctuations peaks at k = as \Uk(t)\ 2 oc 
k~ 2 , leading to long-range correlations, a phenomenon called "dynamical Bose- 
Einstein condensation" in Ref. [7[ and investigated further, for finite volume, in 



Ref. ||18||. We show an example of the momentum distribution in Fig. 5, but we 
have not studied the phenomenon in detail. 

The relation between the asymptotic value as t — > oo for <p(t) and the initial 
amplitude 0o is displayed in Figs. 6 to 8, for T = 1, 2.5 and 3. We compare the 
data with our generalization ( |4.12| ) of the empirical formula fl5~4 ) given in Ref. 



0. The data are obtained by averaging over the second half of the time interval. 
The agreement is excellent, except at the phase boundary where the averaging 
converges slowly. 

As an illustration of the behavior of the system in the unstable region in the 
one-loop approximation we show, in Fig. 9, the evolution of the field amplitude, 
and, in Fig. 10, the exponential behavior of the fluctuation integral and of the 
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effective mass squared A4 2 (t) of the classical field. 

The behavior of the system in region III is displayed in Figs. 11 and 12. The 
amplitude <p(t) is seen to decrease to zero. The decrease is powerlike, not expo- 
nential, a phenomenon called anomalous relaxation in Ref. f6j. Fig. 12 shows 
the squared mass A4 2 (t) which is seen to converge to an asymptotic value M^. 



The sum rule for this asymptotic value, Eq. (|4.12|) , is compared to the data in 



Fig. 13 for T = 1.5, 2.5 and 4. The agreement is again excellent. 

We have not presented the results for the pressure and the ratio of pressure 
and energy which varies between for a nonrelativisticce[Dce[Dce[Dce[Distic and 
1/3 for an ultrarelativistic ensemble. Here these are dominated, already at T = 1, 
by the purely thermal contributions, so that the fluctuations generated by the 
motion of the field <p{t) are relatively unimportant. 



6 Conclusions and Outlook 

The dynamical exploration of the quantum states of the O(N) A$ 4 theory in the 
limit N — > oo has been extended here to finite temperature. We have performed 
numerical simulations with various initial fields 0o = 0(0) and initial masses itlq 
related by the gap equation, and for various temperatures T. Depending on the 
initial conditions we find, in analogy to computations at zero temperature [TJ, 
final states with restored O(N) symmetry and final states for which the symmetry 
is spontaneously broken. The resulting phase diagrams resemble typical phase 
diagrams of thermodynamical systems, with the temperature and an external 
variable as parameters. Instead of, e.g., the magnetic field or the pressure we 
have here the initial value 4>o as external parameter. While the initial states are 
thermal states, the final states are not. 

We have generalized two empirical formulae, the relation between the initial 
and asymptotic field amplitudes in region II, and the formula for the asymptotic 
value of M. 2 {t) in region III to finite temperature, extending the plausibility argu- 
ments given in [[/]. While we have not been able, either, to derive these formulae, 
the way of generalizing them may give some clue for such a derivation. Both 
relations are linked as they give the same formula for the boundary between re- 
gions II and III, though the arguments for their heuristic derivation are seemingly 
different. Furthermore, it is clear that both of them are based on the early time 
behavior. Obviously the fluctuations have to be included up to order V 2 (t) in 
a perturbative expansion. At T = these terms are essentially absorbed into 
renormalization constants, so that the purely classical behavior prevails. One 
may also formulate the modifications at finite temperature in terms of tempera- 
ture dependent masses and couplings. It is the role of the large- N quantum back 
reaction to transmit the early time behavior into the late time one. 

Unfortunately there are many interesting models for which the large-iV ap- 
proach is not possible or not adequate. The one-loop approximation, on the other 
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hand, can be applied in general. However, it shows features that seem to make 
it obsolete for describing nonequilibrium phenomena, especially for theories with 
spontaneous symmetry breaking. As an illustration we have shown the typical 
behavior of a spontaneously broken A$ 4 model in the one-loop approximation. 
The system does not reach a stationary state at late times: the effective mass of 
the classical field diverges exponentially, while the effective mass of the quantum 
fluctuations is and stays negative. This is due to the lack of the quantum back re- 
action onto the fluctuations. The fact that one finds such a pathological behavior 
may, however, indicate the correct physics and is not necessarily a consequence 
of an inadequate approximation. It is known that the system is indeed unstable 
for spatially constant static fields, it is an instability with respect to formation of 
domains |TJ| . For space dependent fields like minimal bubble configurations the 



one-loop approximation to the effective action does not display any unplausible 
features [^, [Tj| |T7| , though the effective potential is complex in the unstable re- 
gion. So it is not clear whether the "taming" of the instability introduced by the 
large- N approximation necessarily improves the understanding of the physics. 

In this situation it is certainly very important to develop new approaches to 
the evolution of quantum systems for theories with spontaneously broken sym- 



metry 19, 20 1 . There are indications in a large- N quantum mechanical system 



[20|] that the large- N limit may be misleading, as the next-to-leading corrections 



become large especially at late times. It is not clear, however, what the impact 
of these results on quantum field theory will be. One of the problems is that, in 
contrast to the large- N and one loop approximations, alternative wave function- 



al pose problems with renormalization [19]. This is not only a technical problem. 



It is connected (trivially) to the fact that the higher the dimension of space, the 
more the ultraviolet behavior of the system will be important. 

We think nevertheless, that a good understanding of the leading order approx- 
imation may improve the understanding of the corrections. That these become 
large at late times is not too surprising, it is therefore even more important to 
realize (once more) that the late-time behavior is related to the early-time be- 
havior, which will therefore set the initial conditions for other approximations as 
well. This should apply to the phase structure as well. 
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A Perturbative expansion 

The mode functions U k (t) with the initial conditions introduced in section 2 
satisfy the integral equation 

oo 

U k (t) = e-^°* + J dt'A ktTet (t - t')V(t')U k (t') , (A.l) 
o 

with 

A fc , rct (t - t') = 0(t - t') sin (u k0 (t - t')) . (A.2) 

We separate U k (t) into the trivial part corresponding to the case V(t) = and 
a function fk{t) which represents the reaction to the potential by making the 
ansatz 

£/*(*) = e- <w *°'[l + /*(*)]. (A.3) 
fk{t) satisfies then the integral equation 

t 

h{t) = J dt'A k , ct (t - t')V(t')[l + f k (t')}e^-^ , (A.4) 
o 

and an equivalent differential equation 

h(t) - 2iu ko f k (t) = -V(t)[l + f k {t)} , (A.5) 

with the initial conditions f k (0) = f k (0) = 0. We expand now f k (t) with respect 
to orders in V(t) by writing 

/*(*) = + (*) + ■■■ (A.6) 

= /^W+ZfW, (A.7) 

where f k (t) is of n'th order in V(t) and f k (t) is the sum over all orders begin- 
ning with the n'th one: 

oo 

/i n) (*) = £/i°(0- (A.8) 

l=n 

The are obtained by iterating the integral equation (|A.4j) or the differential 

equation (|A.5| ). The function f k (t) is identical to the function f k (t) itself which 

is obtained by solving ( |A.5| ). The function f k (t) can again be obtained by 
iteration via 

t 

ff\t) = J dt'A k , ct (t - l!)V$)ft\i!)e**«W . (A.9) 
o 
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The integral equations can be used in order to derive the asymptotic behavior 
as u k0 — > oo and to separate divergent and finite contributions. This has been 
described previously in extenso [O]. We illustrate the procedure by calculating 
the relevant leading terms for fjpit). We have 

t 

fk\t) = ^— I dt'(eM^ko(t - 0) - WO • (A.IO) 
Integrating by parts we obtain 
fi\) = ldt'V(t')--± T -V(t) + — r [dt'exp(2iu k0 (t-t'))V(t') , (A.ll) 

For the expansion of the fluctuation integral T{t) we need the real part of f k 
for which we find 

1 1 r 

Re h£\t) = --r^-V(t) + —J- dt'cos(2uj k o(t-t'))V(t') . (A.12) 



o 



The second term decreases at least as co k Q- In terms of the perturbative expansion 
for the functions f k we can the mode functions appearing in the fluctuation 
integral as 

|t/ fc | 2 = l + 2Re/i TT +|/i TT | 2 . (A.13) 



Using Eq. ( A.12 ) the leading behavior of this expression is 

1 

J k0 



1 + 2Re/i 1} + |if| 2 ^ i _ __ V (t) . (A. 14) 



Similarly the integrand of the energy density and pressure can be expanded [|13[ . 
As these are more divergent, the calculations require more integrations by parts 
in order to single out the leading powers in uj k o and they become more involved. 
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Figure Captions 



Fig. 1: Phase diagram in the (f^ — T plane. 
Fig. 2: Phase diagram in the Wq — T plane. 
Fig. 3: Evolution of classical field in region II. 
Fig. 4: Evolution of M. 2 (t) in region II. 

Fig. 5: The momentum spectrum for T = 1 at t = 75 displaying "dynamical 
Bose-Einstein condensation", with a fit sin 2 (kt)/k 2 . 

Fig. 6: Late time amplitude <p(oo) vs. initial amplitude 0o for T = 1 (asteriscs), 
compared with Eq. ( |4.12j) (solid line). 
Fig. 7: The same as Fig. 3 for T = 2. 
Fig. 8: The same as Fig. 3 for T = 3. 

Fig. 9: Evolution of the classical field in the one-loop approximation. 

Fig. 10: The fluctuation integral (solid line) and M. 2 {t) (dashed line) in the 

one-loop approximation. 

Fig. 11: Evolution of the classical field in region III. 
Fig. 12: Evolution of M 2 (t) in region III. 

Fig. 13: The asymptotic sum rule for A4 2 (t). The data for T = 1.5 (diamonds), 
T = 2.5 (asteriscs) and T = 4 (triangles) are compared to Eq. (|4.19| ) (solid lines). 
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